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ABSTRACT 


Solutions  are  derived  for  the  potential  distributions 
over  one-layer  and  two-layer  random  conductivity  earth  models 
for  the  direct  current  resistivity  method.  The  random 
potential  due  to  a stationary  Gaussian  random  conductivity 
function  is  non-Gaussian  and  non-stationary . The  ensemble 
and  sample  statistics  of  the  random  potential  fields  are 
examined . 

The  potential  field  due  to  dipole  excitation  is  more 
sensitive  to  variations  in  conductivity  in  the  subsurface 
than  that  due  to  monopole  excitation.  Random  variations  in 
conductivity  at  depth  are  difficult  to  detect. 

Representative  curves  of  the  apparent  resistivity  and 
kernel  functions  are  presented  to  show  the  effects  of  a 
random  conductivity  profile  on  them.  These  effects  are 
appreciable  and  result  in  significant  errors  in  the  inter- 
pretation of  resistivity  data.  However,  noise  of  this  kind 
is  difficult  to  distinguish  from  the  signal  itself,  especially 
in  the  presence  of  measurement  noise.  Nevertheless,  variations 
in  the  conductivity  profile  should  be  recognized  as  an 
additional  source  of  error  in  the  interpretation  of  resistivity 
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INTRODUCTION 

One  of  the  principal  objectives  of  electrical  methods 
in  geophysical  prospecting  is  the  study  and  interpretation 
of  resistivity  distributions  in  the  earth.  To  facilitate 
the  interpretation  of  data  simple  earth  models  are  assumed. 

The  measured  potential  field  is  compared  to  the  theoretical 
field  computed  for  the  models.  Interpretation  of  subsurface 
electrical  structure  is  based  on  the  model  which  provides  the 
best  fit.  In  order  to  reduce  the  problem  to  manageable  pro- 
portions, it  has  been  standard  practice  to  consider  only 
discrete  (or  deterministic)  model  parameters.  Such  assumptions 
have  generally  been  adequate  for  data  interpretation.  However, 
it  is  useful  and  instructive  to  utilize  models  which  take 
into  consideration  the  inherent  variability  of  the  earth’s 
physical  properties. 

As  a step  in  this  direction,  this  study  is  directed  at 
the  problem  of  determining  the  potential  distribution  which 
arises  from  direct  current  flow  in  a medium  in  which  the  con- 
ductivity varies  randomly  with  depth.  The  noise  which  is 
consequently  introduced  into  the  system  is  best  handled 
statistically . 

An  earth  model  with  a random  conductivity  profile  has 
been  previously  considered  in  only  one  paper  (Naidu,  1970) 
in  the  readily  available  geophysical  literature.  In  this 
paper  only  the  one-layer  isotropic  case  was  treated  in  any 
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detail.  This  study  extends  the  scope  of  the  original  paper 
by  examining  in  considerable  depth  various  aspects  of  the 
statistical  properties  of  random  earth  potentials.  An 
analytic  solution  to  the  more  useful  and  practical  two-layer 
case  is  also  presented.  The  effect  of  anisotropy  on  the 
random  conductivity  model  is  examined.  The  implications  of 
the  random  distribution  for  resistivity  inversion  are  also 
considered. 

The  primary  purpose  of  this  dissertation  is  to  study  in 
some  detail  the  direct  (or  forward)  problem  without  an  under- 
standing of  which  all  the  information  contained  in  field  data 
cannot  be  effectively  utilized.  It  is  of  considerable 
practical  importance  to  be  able  to  judge  quantitatively  the 
effects  of  the  random  properties  of  the  conductivity  distri- 
bution in  the  earth  on  the  direct  current  resistivity  method. 
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THE  RANDOM  CONDUCTIVITY  PROBLEM 


Direct  current  resistivity  methous  used  in  geophysical 
exploration  consist  essentially  of  the  measurement  and  inter- 
pretation of  electric  potentials  on  the  surface  of  the  earth. 
These  potentials  result  from  the  introduction  of  an  electric 
current  into  the  ground  through  various  electrode  arrays. 

The  observed  potential  field  is  a function  of  the  resistivity 
distribution  in  the  subsurface.  The  interpretation  of  the 
acquired  data  in  terms  of  subsurface  electrical  structure  is 
an  inverse  (or  inversion)  problem.  For  this  inversion  tc  be 
possible,  the  solution  to  the  direct  (or  forward)  problem  is 
required,  that  is,  given  a model  of  the  earth,  one  needs  to 
compute  the  corresponding  potential  distribution  on  the  surface. 

Statement  of  the  Problem 

The  simplest  and  most  widely  used  model  assumes  a half- 
space composed  of  discrete  horizontal  layers  with  sharp 
discontinuities  in  the  electrical  properties  of  the  medium 
at  the  layer  boundaries.  Each  layer  is  homogeneous  and 
assumed  to  be  either  isotropic  or  anisotropic.  Isotropy  is 
t.he  more  common  assumption  however.  The  solution  to  this 
basic  resistivity  problem  is  well  documented  in  the  literature, 
for  example  Grant  and  West  (1965),  Keller  and  Frischknecht 
(1966),  Van  Nostrand  and  Cook  (1966)  and  Bhattacharya  and 
Patra  (1968).  Solutions  have  also  been  published  for  models 
in  which  the  conductivity,  instead  of  being  constant  over 
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discrete  vertical  intervals,  is  a continuous  function  of  depth 
(Slichter,  1933,  Langer,  1933,  Meinardus,  1967  among  others). 

In  this  study  the  problem  is  generalized  by  assuming 
that  the  conductivity  is  not  a deterministic  function  of 
depth  z but  is  instead,  a random  function  of  z which  is, 
nevertheless,  closely  clustered  around  its  deterministic  value. 
It  was  noted  by  Keller  (1968)  that  resistivity  probability 
density  curves  usually  represent  multimodal  log-normal  dis- 
tributions; this  is  due  to  the  presence  of  several  lithologic 
types.  Tc  make  the  problem  tractable  mathematically,  only 
unimodal  normal  distributions  of  conductivities  are  considered 
here . 

Consider  a model  (Figure  1)  in  which  the  first  layer  (the 
overburden)  of  thickness  h has  a random  conductivity  profile 
given  by 

^(z)  - o1°(l  + en(z)  ) (2-  1) 

where  a^,  a constant,  is  the  deterministic  (in  this  case,  the 
mean)  conductivity,  e is  a small  parameter  (s  <<  1)  and  n(z) 
is  a zero-mean  Gaussian  random  function  with  a prescribed  co- 
variance  function  C.  Thus, 

< n(z)  > - 0 (2-  2) 

and  < n(z1)  n(z2)  > =•  C( | zx  - z2|)  (2-  3) 

where  < • > denotes  ensemble  average  or  expected  value. 

The  second  layer  (the  basement)  is  assumed  to  have  a 
constant  conductivity  a2<  Both  layers  are  assumed  to  be 
isotropic.  The  first  layer  is  isotropic  only  on  the  micro- 


Figure  1.  The  two-layer  random  conductivity  model. 


T-1950 


6 


scopic  scale;  macro-anisotropy  is  implicit  in  a random  con- 
ductivity profile. 

The  problem  is  to  determine  the  potential  distribution 
which  would  result  from  the  introduction  of  a direct  current 
of  strength  I into  the  medium  described  above  through  an 
electrode  located  on  the  surface.  The  solution  is  derived 
by  solving  the  applicable  differential  equation  with  the 
appropriate  boundary  conditions. 


Boundary  Conditions 

Laplace's  equation,  in  cylindrical  coordinates,  for 
inhomogeneous  isotropic  media  is 

_ / _ v 3 2 V ,13V,  da(z)  3V  , _,_,3ZV  _ x 

0<Z)77  r 'Sr  * <’(Z)77  0 ■ ( 4) 

3r  32 

This  equation  applies  to  both  layers  and  V is  the  electric 

potential  in  each.  Let  Va(r,z)  and  Vb(r,z)  be  the  potentials 

in  the  first  and  second  layers  respectively. 

The  boundary  conditions  to  be  satisfied  can  be  summarized 


as 

follows : 

(a) 

(r,z)  - 0 

as 

r 

•+•  00  , i * a , b 

(2-  5) 

(b) 

Va  - 

I 

as 

r 

*►  0 at  z ■ 0 

(2-  6) 

a 

2tt0j  (o)r 

(c) 

vb  - o 

as 

z 

-*■  co 

to 

i 

-o 

3Va 

f o 

I 

r f 0 

(d) 

a 

sa 

< 

(2-8) 

3z 

z » 0 

r " 

i 

r ■ 0 

(e) 

Va  “ V.  at 

a b 

2 ■ 

h 

(2-  9) 
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(f) 


ax(h) 


3Va 

TT 


z = h 


3Vt 

<Vh>  Tz 


(2-10) 


Properties  of  the  Random  Function 

Before  proceeding  with  the  solution,  it  is  pertinent  to 
review  the  properties  of  the  stationary  Gaussian  random 
function  n(z)  which  will  be  useful  in  formulating  the  solution 
to  the  boundary  value  problem. 

The  integral  canonical  representation  (Pugachev,  1965, 
p.  309)  of  n(z)  is 

n(z)  = / N(w)ela)Z  dw  (2-11) 


where  N ( a,)  is  white  noise  of  the  transform  variable  u,  the 
intensity  of  which  is  equal  to  the  spectral  density  S(uj  ) of 
the  random  function  n(z).  In  other  words,  the  random  process 
N(to)  is  zero-mean  and  uncorrelated  (see  also  Koopmans,  1974), 
that  is, 


< N(u1)N(uj2)  > - S(ux)  6(w1  - u2)  (2-12) 

where  <5  is  the  Dirac  function. 

It  is  further  noted  that 

N(  a) ) - -i-  / n(z)e"iuz  dz  . (2-13) 

2t  -» 
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If  (z^-Zg)  is  denoted  by  x , it  can  be  shown  (Pugachev,  1965) 


that 


and 


C(x) 


/ S(w)e 


lux 


do) 


(2-14) 


s<">  - £ 7 cuxr1"* 


dx  (2-15) 

where  C(x)  is  the  covariance  of  n(z).  Thus  C(x)  and  S(<u) 


form  a Fourier  transform  pair. 

Boundary  condition  (2-8)  requires  that  n be  an  even 
function  of  z.  This,  in  turn,  implies  that  N(w)  is  an  even 
function  of  uj  . These  symmetry  relations  will  simplify  con- 
siderably the  equations  above  so  that  only  real  terms  remain 
as  the  imaginary  ones  drop  out. 

Assume  a covariance  function  of  the  form 

C(x)  » q2exp{-x2/(4a2)  ) (2-16) 

o 

where  q is  the  variance  of  n(z)  and 

a is  an  arbitrary  constant  which  determines  the  degree  of 

correlation.  The  exponential  term  in  equation  (2-16)  is  the 

correlation  function  of  n(z).  Figure  2 shows  the  function 

C(x)  for  various  values  of  the  parameter  cc. 

Since  the  covariance  function  C of  a real  stationary 

random  function  is  an  even  function  of  x,  its  spectral  density 

S is  also  an  even  function.  Corresponding  co  C(x)  as  given 

above  we  have  the  spectral  density  function 

„ 2 2 
a 2 -a  a) 

■ — q e 

/tt 


S(w) 


(2-17) 
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General  Solution 

To  solve  the  two-layer  stochastic  problem,  as  with  the 
deterministic,  the  variables  are  first  separated: 

V(r,  z)  = R(r)P(z) 

Equation  (2-4)  then  becomes 

R"(r ) + ^R'(r)  + \2R(r)  * 0 (2-18) 

and 

(o(z)  P'(z)  )’  - X2o(z)  P(z)  = 0 (2-19) 

where  X is  the  separation  constant  and  primes  denote  differen- 
tiation with  respect  to  the  parameters. 

The  equation  in  r is  Bessel's  equation  and  its  solution 

is  JQ( Xr)  for  both  layers.  The  solution  to  equation  (2-19) 

“ X z 

for  the  second  layer  is  e . For  the  stochastic  first  layer, 
however,  the  solution  to  (2-19)  is  more  mathematically  involved 
and  the  reader  is  referred  to  Appendix  A for  the  complete 
development  of  the  solution.  Only  the  final  expressions  will 
be  presented  here. 

The  potential  on  the  surface,  to  a first-order  approximation, 
can  be  expressed  as 

V(  r , o)  = Vo(r,  o)  + eV^r,  o)  + 0(e2)  (2-20) 

where  r is  the  distance  of  the  point  of  measurement  from  the 
current  source.  The  first  term  Vq  is  recognized  as  the  potential 
due  to  a "deterministic"  two-layer  case  in  which  both  layers 
have  constant  conductivities.  Thus 
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V ( r , o ) = ~~  / K (A)J  (Xr)dA 


a^°  o 


(2-21) 


where  A ^ i— 

4" 


Kq(A) 


1 + ke 
1 - ke 


-2Ah 

^2lh 


(2-22) 


The  reflection  coefficient  k is  defined  as 

o o 

, A °1  - a2 

7"° — 7^ 

°1  + a2 

For  convenience,  the  function  Kq  will  be  referred  to  as  the 
deterministic  kernel  (corresponding  to  the  so-called  Slichter 
kernel  in  deterministic  models)  and  VQwill  be  called  the 
deterministic  potential. 

The  term  can  be  viewed  as  a first-order  perturbation 

of  V due  to  the  fact  that  the  conductivity  of  the  first  layer 
o 

is  not  a constant,  but  is  a random  function  of  depth. 

In  Appendix  A,  the  following  solution  is  derived: 


Vl( r , o)  =■  ~ / Kx(l)  J0(Ar)dX 

o j o 


(2-23) 


2 00 

where  K.(A)  = - — — — / — du 

1 dou 


, 8 X2(  l-k2)e"2Xh  " N cos  u)h 

+■  5 j uu> 

ip  o u 


4A( l+k2)e"2Ah  * N w sin  uh 

1 1 J 1 ■ ' ' ~ ■*" 


0 


(2-24) 


L...C ■ 'Ml i*i ifariUhii  ■ . a-  -■‘J  , - L.. 


T-1950 


12 


For  brevity  m notation  the  following  symbols  have  been 
used : 


1 - ke 


-2Xh 


v 

y 


-2Xh 


1 + ke 
w2  + 4X2 


It  is  convenient  to  refer  to  and  as  the  stochastic 
(or  random)  potential  and  the  stochastic  (or  random) 
kernel  respectively. 

Equations  (2-21)  through  (2-24)  together  represent  the 
general  solution  to  the  two-layer  stochastic  boundary  value 
problem. 


Reduction  to  the  One-Layer  Case 

Let  the  stochastic  first  layer  extend  to  infinity  in 
the  z direction.  The  model  then  becomes  a half-space  whose 
conductivity  varies  randomly  about  its  mean  (or  equivalently, 
a randomly  layered  semi-infinite  medium). 

The  solution  to  this  boundary  value  problem  is  readily 
obtained  from  the  general  solution  above  by  letting  the  over- 
burden thickness  h tend  to  infinity.  Thus 


K0  (X) 
VQ(r,  o) 


K:(a) 


1 


5 * N 

-8X^  / y du) 

o 


(2-25) 

(2-26) 


(2-27) 
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In  this  special  case  KQ  and  VQ  are,  respectively,  the  kernel 
and  potential  functions  of  a homogeneous  isotropic  half-space. 

is  the  perturbation  potential  which  arises  from  the  random 
nature  of  the  conductivity  profile. 

The  equations  above  are  in  agreement  with  those  derived 
by  Naidu  (1970)  for  the  one-layer  case.  The  anisotropic  one- 
layer  stochastic  problem  is  solved  directly  in  Appendix  B 
using  a procedure  different  from  Naidu's.  In  addition,  the 
concept  of  anisotropy  is  carried  further  here  than  in  the 
original  paper. 
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THE  ONE-LAYER  CASE 


Considerable  insight  into  the  random  conductivity 
problem  may  be  gained  by  examining  first  the  statistical 
properties  of  the  single  layer  case.  To  simplify  the  dis- 
cussion, it  will  be  confined  mainly  to  the  isotropic  case. 
Anisotropy  is  dealt  with  separately  at  the  end  of  the  chapter. 


Monopole  and  Dipole  Potentials 

All  the  expressions  for  the  potentials  which  were  pre- 
sented in  the  previous  chapter  apply  to  a monopole  source. 

If,  instead,  the  source  is  a dipole  (Figure  3),  it  can  be 
shown  that  we  merely  have  to  replace  A by  B and  JQ(Xr)  by 
XJ-^Xr)  in  all  previous  equations.  The  kernel  functions,  being 
independent  of  electrode  configurations,  remain  unchanged. 

Thus,  denoting  dipole  potentials  by  U,  to  distinguish  them 
from  monopole  potentials  V,  we  have 


UQ(r,  o) 


B 


2 o 
r ax 


U x ( r , o) 


~o“  7 K,  ( X ) XJ.(Xr)  dX 

al  o 1 1 


where  is  given  in  equation  (2-27)  and 


B 


^ II  cos  8 


(3-1) 

(3-2) 


£ = half-length  of  dipole  source 


0 = azir  . of  the  point  of  measurement  from  the 


dipole  . 
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For  the  sake  of  brevity  in  notation  all  explicit 
references  to  the  parameters  of  the  potentials  will  be 
suppressed;  unless  otherwise  stated,  they  are  understood 
to  be  (r,  o) . 

After  normalizing  the  stochastic  potentials  by  their 
respective  deterministic  potentials,  reversing  the  order  of 
integration  and  integrating  once  we  obtain 

Vx  - — - 2 / { | ^ Q0(iy)  - 1 > N(w)  dw 

Vo  0 


= 2 / g^(wr)  N(w)  dw 

o 


(3-  3) 


and 


o t— 

00 

= 2 / g2(wr)  N(w)  dw 


dw 


“i  ■ Sr1  " 2 * ['I1'2  { 1 - 5 vr5  1 - 0 N(u) 

(3-  4) 


A 0 

where  Q^(x)  = I^(x)  - L^(x) 


I (x)  = modified  Bessel  function  of  order  i 


L^(x)  = modified  Struve  function  of  order  i 


For  details  of  the  derivation  of  the  equations  above 
the  reader  is  referred  to  Appendix  B. 

Ensemble  Statistics 

Evidently  the  normalized  random  potentials  ana 
are  non-Gaussian  random  functions.  It  is  instructive  to 
examine  their  ensemble  statistics. 
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Their  means  are  obviously  zero  because  the  random  function 


N(oj)  is  zero-mean.  The  covariance  of  is 


00  oo 


<Vl^rl^Vl^r2 }>  * f g1(u)1r1)g1(iD2r2)  <N(u1)N(uj2)>  du-^do^ 


0 0 

00 


*1 v 1 ly&lv22' 


4 / g1(ujr1)  g^wrg)  S(w)  du> 

, 2 « 2 2 
— • ■*-  f g1(<*>r1)  g1(uar2)e_  ^ 


doj 


(3-  5) 


It  is  to  be  noted  that  equations  (2-12)  and  (2-17)  have  been 
used  in  the  derivation  above. 

The  normalized  covariance  is  defined  as 


Rl(rl*  r2>  = 


< Vx(ri)  V1(r2)  > 


{<Vl(rl)  Vl(rl}>  <Vl(r2)  Vl(r2}>  ]i  (3'  6) 

For  dipole  excitation  the  corresponding  expressions 


are 


2 » 2 2 
<ul(ri)ui(r2)>  = f g2(UJr1)g2(aJr2)e‘a  w 

^ vTT  0 *“ 


da) 


(3-  7) 


and 


R2^rl’  r2^ 


< U1(r1)  Ux(r2)  > 


(<ul<rl>  U1(r1)>  <ul(r2)  t^CrpM*  ( 3.  g) 


The  Q1  functions  were  evaluated  by  means  of  polynomial 
approximations  and  power  series  expansions  (Abramowitz  and 
Stegun,  1965,  p.  378  and  p.  498).  As  the  argument  increases 
from  zero  Qq  decreases  raonotonically  from  1 . 0 to  zero  while 
Q-^  increases  raonotonically  from  zero  to  its  asymptotic  value  of 
2/tt  . The  exponential  term  in  equations  (3-5)  and  (3-7) 
ensure  rapid  convergence  of  the  integrals.  Using  numerical 
integration  techniques  the  normalized  covariance  functions 
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and  were  computed  for  five  representative  values  of  a and 
for  various  distances  r from  the  source. 

Some  typical  results  are  shown  in  Figures  4 through  7. 

From  equation  (2-16)  and  Figure  2,  it  is  obvious  that  a 
smaller  value  of  the  parameter  a implies  that  the  random 
function  n(z)  is  less  correlated  and  vice  versa.  Figures  4 
and  5 show  that  as  a decreases  the  covariance  curves  become 
narrower,  that  is,  the  normalized  random  potentials  become 
less  correlated.  Hence,  it  can  be  concluded  that  as  the 
correlation  of  the  random  conductivity  function  decreases,  the 
correlation  of  the  normalized  random  monopole  potential  decreases 
too. 

Similar  conclusions  can  be  drawn  for  the  normalized 
random  potentials  which  result  from  dipole  excitation 
(Figures  6 and  7). 

Sample  Statistics 

The  study  of  the  sample  statistics  of  the  random 
potentials  requires  the  generation  of  a large  number  of  sample 
realizations  of  the  random  potentials  for  both  the  monopole 
and  dipole  cases.  (A  sample  realization  is  an  observation  on 
a random  process.)  In  order  to  do  this  the  random  potentials 
have  to  be  expressed  in  terms  of  the  function  n(z)  explicitly. 

For  the  raonopole  case,  it  is  recalled  from  equations 
(2-23)  and  (2-27)  that 

V = A*  7 - 8A2J  (Ar)  / 


dcu  dA 


(3-  9) 


Figure  4.  Normalized  covariance  of  monopole  random  potential 
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♦ o 

where  A = A /a^ 

2 2 

and  u = w + 4A 

Substituting  for  N(oj)  with  equation  (2-13)  and  changing  the 

order  of  integration  we  have 

$ 00  00  00 

V,  - - — / \2J  (Ar)  / n(z)  f — dtadzdA  . (3-10) 

1 n o 0 o oh 

Evaluating  the  third  integral  analytically  (Dwight,  1961) 
we  obtain 

00  00  n \ 

V,  = -2A*  I n(z)  f Xe~2Xz  J ( Ar)  dA  dz  (3-11) 

1 o o ° 

Differentiating  the  Lipschitz  integral  (Watson,  1966) 

/ e“2Az  J ( Ar ) dA  = (r2  + 4z2)  * 

c ° 

with  respect  to  z the  second  integral  in  equation  (3-11) 
can  also  be  evaluated  analytically.  Thus, 


v - -4A*  t ’l<|>  Z df  ,/, 

1 o (r*  + 4 z“) 

= -4A*  M(r) 

When  the  dipole  case  is  similarly  handled  we  have 

U „ _12B*  r 1 
1 o (r2  + 4z2)5/2 

= -12B*  D(r) 


(3-12) 


(3-13) 


where 


T-1950 


24 


The  stationary  Gaussian  random  function  n(z)  was 
generated  by  standard  procedures  (see,  for  example,  Newman 
and  Odell,  1971  and  Hemmerle,  1967).  Essentially,  Gaussian 
weights  were  applied  to  unit  normal  random  variables  obtained 
via  central  limit  convergence  from  uniform  random  variables. 
The  latter  were  supplied  by  a random  number  generator.  The 
subroutine  for  generating  n(z)  was  modified  from  a version 
prepared  by  Barakat  (personal  communication,  1976). 

It  is  evident  that  the  denominators  in  the  integrals 
M(r)  and  D(r)  will  ensure  rapid  convergence.  Two  hundred 
sample  realizations  for  each  of  several  selected  values  of  r 
and  ct  were  computed  for  both  functions  M(  r)  and  D(r).  The 
results  were  tabulated  in  the  form  of  histograms,  examples 
of  which  are  presented  in  Figures  8 and  9 for  the  mcnopole 
and  dipole  cases  respectively. 

The  computed  means  are  generally  less  than  0.01.  These 
are  close  enough  to  the  theoretical  value  of  zero  within  the 
limits  of  sampling  error.  The  standard  deviations  are  of  the 
order  of  0.1  and  the  coefficients  of  skewness  (normalized 
third  product  moments)  are  not  greater  than  about  0.1  in 
magnitude . 

All  distributions  of  the  random  monopole  and  dipole 
potentials  generated  for  various  values  of  r and  a were  sub- 
jected to  the  chi-square  goodness  of  fit  test  (Walpole  and 
Meyers,  1972,  for  example).  It  was  found  that  at  the  5% 
level  of  significance  the  Gaussian  distribution  provides  a 


0.15 


monopole  random  potential 
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good  fit  for  all  the  distributions  of  random  potentials  that 
were  generated.  The  sample  statistics  therefore  indicate 
that  the  normalized  random  potentials  and  are  normally 
distributed,  at  the  5%  level  of  significance,  about  their 
mean  (deterministic)  values. 

Sample  Realizations 

It  is  informative  to  see  how  the  random  function  n(z) 
affects  the  measured  voltage,  the  kernel  function  and  the 
apparent  resistivity.  Some  examples  are  shown  here. 

As  noted  earlier,  the  total  potential  for  monopole 
excitation,  to  a first-order  approximation,  is 

V(r)  = VQ(r)  + e V^r) 

* - e 4A*  M(r) 

Normalizing  by  A*  the  expression  above  becomes 

V*(r)  - i - 4 e M ( r ) . (3-14) 

Similarly,  for  the  dipole  case,  after  normalizing  by  B* 
the  equation  is 

U*(r)  =>  -9  - 12eD(r)  . (3-15) 

Four  representative  sample  realizations  of  the  total 
potential  measured  over  a random  medium  for  both  source  arrays 
are  shown  in  Figures  10,  12,  14  and  16.  In  the  figures,  the 
solid  lines  represent  the  ''normal”  potential  measured  over  a 
homogeneous  half-space.  The  curves  with  circles  superimposed 
show  what  the  measured  voltages  will  be  if  the  random  con- 
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Figure  13.  Sample  realization  No.  2 of  kernel  and  apparent 
resistivity  functions  over  random  conductivity 
model. 
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Figure  15.  Sample  realization  No.  3 of  kernel  and  apparent 
resistivity  functions  over  random  conductivity 
model . 
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Sample  realization  No.  4 of  kernel  and  apparent 
resistivity  functions  over  random  conductivity 
model . 
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ductivity  functions  are  as  shown  at  the  top  of  the  figures. 

The  parameter  a is  arbitrarily  set  at  0.15  (corresponding 
to  a weakly  correlated  n(z)  ).  To  make  the  variations  in  the 
random  voltage  large  enough  to  be  discerned  on  the  graphs,  the 
parameter  e has  been  set  at  a relatively  large  value  of  0.1. 

It  is  apparent  that  in  the  transformation  from  a random  con- 
ductivity distribution  to  the  measured  potential,  high  frequency 
information  has  been  lost.  The  random  potentials  corresponding 
to  each  conductivity  distribution  are  smooth  functions  of  r . 
Herein  lies  the  difficulty  of  distinguishing  noise  in  the 
system  from  the  signal  itself. 

It  is  more  useful  to  examine  the  kernel  and  apparent 
resistivity  functions.  The  stochastic  kernel  given  in 
equation  (2-27)  can  be  expressed  in  terms  of  n(z)  explicitly 
as 

K1(X)  = - 2X£n(z)e'2Az  dz  . (3-16) 

The  "total"  kernel  is 

K( X ) = KQ(X)  + cKx( X) 

* 1 - 2£Xjn(z)e“2Az  dz  (3-17) 

Let  us  examine  the  asymptotes  of  the  kernel  function.  Clearly, 

K1(o)  = 0 

It  is  less  obvious,  but  nevertheless  readily  proved,  that 
K^(  X)  -+■  -n(o)  as  X -*■  » 
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It  should  be  noted  that  in  computing  K^(X,  by  numerical 
integration,  we  need  a mesh  size  Az  that  becomes  progressively 
smaller  as  X gets  correspondingly  larger.  This  is  necessary 
for  a reliable  approximation  to  the  value  of  the  integral. 

The  total  kernel  functions  K(X)  for  the  same  four  sample 
realizations  that  were  used  in  the  previous  figures  are  shown 
in  Figures  11,  13,  15  and  17.  The  deviations  from  the  normal 
kernel,  KQ(X)  = l,  depend  on  the  nature  of  the  random  con- 
ductivity function.  However,  the  most  important  factor  is 
its  value  at  the  surface  n(o).  For  all  practical  purposes  this 
asymptote  is  reached  for  X values  of  about  100. 

The  normalized  apparent  resistivity  measured  with  a 


pole-pole  (or  single  pole)  array  may  be  expressed  as 
*.  . pa^ 


P (r)  * 


= 1 - 4erM(r) 


(3-18) 


where  pa(r)  is  the  apparent  resistivity 

ai 

and  P°  is  the  reciprocal  of  ai°  . 

The  corresponding  expression  for  the  pole-dipole  array  is 


P*(r) 


Pa(r) 


« 1 - 12er  D(r) 


(3-19) 


These  normalized  apparent  resistivities  are  shown  in 
Figures  11,  13,  15  and  17  in  which  they  are  referred  to  as 
"monopole"  and  "dipole"  apparent  resistivity  curves.  For  a 
homogeneous  half-space  these  normalized  apparent  resistivities 
have  constant  values  of  one.  The  deviations  from  this  mean 
value  reflect  the  random  nature  of  the  conductivity  profile. 
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Clearly,  these  deviations  approach  zero  as  distance  r from 
the  source  increases. 

The  corresponding  Wenner  and  Schluraberger  apparent 

resistivities  will  now  be  briefly  examined.  The  normalized 

Wenner  apparent  resistivity  is 
* P (a)  99 

Pw(a)  - — - 2a  / K(X)  [jQ(Xa)-J0(2Xa)]  dX 

p o 

« 1 - 8ea  (M(a)  - M(2a)  ) (3-20) 

where  a is  the  electrode  spacing  (Figure  18).  It  can  be 

r 

£ 

I shown  that 

Pw*  - 1 - en(0)  as  a - 0 

and 

*-*•1  as  a ® 

pw 

The  normalized  Schlumberger  apparent  resistivity  is 

* P ( a ) _ ® 

Ps  (a)  - -— ■? * a'2  / K(X)XJ1(Xa)  dX 

p o 

- 1 - 12ea2D(a)  (3-21 

where  the  electrode  spacing  a is  as  shown  in  Figure  18.  It 
is  obvious  that 

# 

pg  ■*  1 - en(0)  as  a -*•  0 

and 

p * -*•  1 as  a ■*  00 


Schlumberger  Array 


Figure  18.  Wenner  and  Schlumberger  electrode  arrays 
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Thus,  for  small  electrode  spacings  the  resistivity 
measured  with  either  array  will  be  equal  to  the  resistivity 
of  the  near-surface  zone  and  for  large  electrode  spacings 
the  effect  of  the  random  profile  will  be  negligible  since  all 
the  variations  tend  to  be  averaged  out  when  a greater  depth 
of  the  section  is  probed.  These  apparent  resistivity  curves 
would  be  similar  in  shape  to  the  kernel  curves  shown  if  the 
independent  variable  on  the  latter  is  regarded  as  1/a  instead 
of  \ . These  resistivity  curves  resemble  the  standard  two- 
or  three-layer  curves  for  deterministic  models.  Although  the 
deviations  from  the  corresponding  curves  for  a homogeneous 
earth  are  relatively  small,  these  shifts  in  the  curves  for 
small  values  of  a are  significant  enough  to  affect  the 
accuracy  of  data  interpretation.  It  is  more  convenient  to 
examine  this  in  the  next  chapter  when  the  two-layer  model  is 
considered . 

Sensitivity  to  Resistivity  Variations 

An  indication  of  the  relative  sensitivity  of  direct  current 
methods  to  variations  in  the  conductivity  profile  may  be 
obtained  by  studying  the  variance  of  the  random  potentials. 

For  monopole  excitation  the  variance  is 

2 2 2 

<Vx(r)  V1(r)>  = / {g1(wr)  }2e-a  u du> 

/tt  0 


(3-22) 
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For  dipole  excitation,  we  have 

_ 2 oo  2 2 

<U1(r)  Ux(r)>  = / {g2(ur)}2e"a  u dm  (3-23) 

/if  o 

The  variance  at  r»0  for  both  source  configurations  is  equal 

o 

to  2 q as  can  be  readily  verified.  Hence,  the  variance  of 
the  normalized  random  potential  for  both  cases  is  twice  the 
variance  of  the  random  conductivity  function  n . This  is  also 
the  maximum  value  of  the  variance. 

The  variances  of  the  monopole  and  dipole  random  potentials, 

2 

after  normalizing  by  2q  , are  shown  in  Figures  19  and  20 
for  a =»  0.15  and  ot  = 1.5  respectively.  It  is  apparent  that 
the  variances  of  the  random  potentials  decrease  as  measure- 
ments are  taken  farther  from  the  source.  This  indicates,  not 
unexpectedly,  that  the  sensitivity  of  the  measurements  to 
variations  in  the  conductivity  profile  decreases  markedly  with 
distance  from  the  source.  The  variance  is  larger  for  a more 
correlated  random  conductivity  function  (that  is,  when  ct  is 
large).  For  a completely  uncorrelated  function  n(z)  the 
variance  is  zero  except  in  the  vicinity  of  the  source.  In 
such  a case,  the  conductivity  variations  cannot,  theoretically, 
be  detected. 

It  is  interesting  to  note  that  the  variance  of  the  normal- 
ized random  dipole  potential  is  always  greater  than  that  due 
to  a monopole.  This  leads  to  the  conclusion  that  measurements 
made  with  a dipole  source  are  more  sensitive  to  conductivity 
variations  in  the  subsurface  than  those  made  with  a monopole 
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source.  Consequently,  the  signal  to  noise  ratio  for  dipole 
excitation  decreases  more  rapidly  with  distance  from  the 
source  than  that  for  monopole  excitation. 

Anisotropy 

The  resistivity  of  a rock  may  depend  on  the  direction 
in  which  a current  flows  through  it.  One  cause  of  this  is 
the  microscopic  structure  of  the  rock;  the  alignment  of 
mineral  grains  may  permit  current  to  flow  more  readily  along 
the  direction  of  the  bedding  plane  or  schistosity  than  perpen- 
dicular to  it.  This  is  referred  to  as  micro-anisotropy. 

A sequence  of  parallel  layers  each  with  its  own  thickness  and 
resistivity  will  also  give  rise  to  a preferred  direction  of 
current  flow.  In  this  case,  we  refer  to  structural  anisotropy 
or  macro-anisotropy.  For  excellent  discussions  on  the  nature 
and  effects  of  anisotropy  the  reader  is  referred  to  Schlumberge 
et  al . , 1934,  Kunetz,  1966  and  Keller,  1968. 

Let  us  first  examine  the  effects  of  micro-anisotropy. 
Consider  a one-layer  medium  in  which  the  radial  and  vertical 
conductivity  functions  are  given  by 

ar(z)  = ar°  (1  + enr(z)  ) (3-24) 

az(z)  = a°  (1  + enz(z)  ) (3-25) 

in  which  the  constants  a 0 and  v _°  are  the  mean  conductivities 
in  the  r and  z directions  respectively  and  nr  and  n2  are 
zero-mean  stationary  Gaussian  random  functions  of  z with  pre- 


scribed covariance  functions. 
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The  coefficient  of  anisotropy  S is  defined  as 

rz  o 

3 = =^-  • (3-26) 

*a 

z 

Case  I:  If  it  is  assumed  that  a_( z)/a_(z)  is  constant 
_ r Z 

for  all  z , then 

nr(z)  - n2(z) 

In  the  interest  of  brevity  in  notation  the  subscripts  (which 
are  now  superfluous)  for  the  function  n and  its  transform  N 
will  be  dropped.  As  shown  in  Appendix  B,  t ie  normalized  random 
potentials  for  monopole  and  dipole  sources  are 

Vr'  * 2 l - l!  N<“>  d“  (3-27) 

UjCr,  o)  - 2 7 [(|£>2  (1  - f QiC^§>>  - l] »(<■>>  • (3-28) 

These  equations  differ  from  (3-3)  and  (3-4)  for  the  isotropic 
case  only  in  the  fact  that  there  is  now  an  additional  factor 
S . 

The  normalized  covariances  for  various  0 and  r^  values 
are  shown  in  Figures  21  through  23.  Inspection  of  these  curves 
leads  to  the  following  conclusion.  The  larger  the  coefficient 
of  anisotropy  the  more  correlated  is  the  random  potential. 

It  can  also  be  shown  that  the  variances  of  the  random  potentials 
for  a fixed  distance  r increases  with  g. 

It  is  apparent  that  the  term  3 in  the  equations  above 
are,  essentially,  no  more  than  a scaling  factor  on  the  distance 
r.  Hence,  not  much  additional  information  is  gained  by  con- 
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sidering  anisotropy  separately,  as  is  well  known. 

Case  II:  Instead  of  letting  ar  be  random,  let  us 
assume  it  is  constant, 


This  implies  nr(z)  *»  0.  The  monopole  and  dipole  potentials, 
as  shown  in  Appendix  B,  are  half  the  corresponding  values 
for  Case  I.  It  is  therefore  unnecessary  to  pursue  this 
case  any  further . 

Macro-anisotropy : Before  leaving  the  subject  of 
anisotropy,  it  is  pertinent  to  note  that  even  in  the  so- 
called  "isotropic"  case,  macro-anisotropy  is  implicitly 
present  when  a random  conductivity  profile  is  assumed. 

This  is  illustrated  in  the  following  manner.  The  average 
transverse  resistivity  of  a succession  of  beds,  in  the 
limit,  as  bed  thicknesses  tend  to  zero,  can  be  expressed 
in  the  form 


pt  - 

where  p(z)  is  the  resistivity  as  a function  of  depth  and 
H is  the  total  thickness  of  the  sequence. 

On  the  other  hand,  the  average  longitudinal  conductivity 
of  the  sequence  is  given  by 

, zjdz 
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where  o(z)  is  the  conductivity  as  a function  of  depth. 

The  coefficient  of  macro-anisotropy  is 

b - ■'0^7 

This  value  is  extremely  close  to  1.0  for  the  conductivity 
distributions  that  have  been  used  here.  Even  for  the 
relatively  large  value  of  0.1  for  the  parameter  e , this 
coefficient  of  anisotropy  is  only  about  1.005  on  the  average. 
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THE  TWO- LAYER  CASE 


The  randomly  conductive  overburden  problem  proposed  in 
this  study  arises  in  nature  when  a sedimentary  basin  is  under- 
lain by  a basement.  One  can  conceive  of  a sequence  of  thinly 
layered  sediments  with  slightly  different  electrical  character- 
istics for  each  layer.  The  conductivities,  however,  deviate 
only  marginally  from  the  mean  of  the  sedimentary  section  as 
a whole.  Frequently  the  "electrical"  basement  will  have  a 
much  higher  resistivity  than  the  overburden  (that  is,  the 
reflection  coefficient  k is  positive) . For  all  practical 
purposes  the  basement  resistivity  is  constant  relative  to  the 
variability  in  the  more  conductive  overburden.  However,  one 
can  conceivably  have  a situation  in  the  field  in  which  the 
basement  is  more  conductive.  The  question  then  naturally 
arises  as  to  whether  the  model,  as  it  is  assumed,  is  applicable. 
It  would  seem  to  be  more  appropriate  in  this  case  to  have  the 
variability  in  the  more  conductive  basement.  It  should  be 
pointed  out,  however,  that  shallower  layers  have  a greater  effect 
on  the  surface  measurements  than  deeper  layers. 

In  any  case,  the  statistical  properties  of  the  random 
overburden  model  for  both  positive  and  negative  reflection  co- 
efficients will  be  examined,  bearing  in  mind  that  the  case  of 
a more  resistive  basement  is  generally  of  greater  interest 
from  a practical  viewpoint.  It  is  appropriate  to  stress  here 
that  although  this  is  referred  to  as  a "two-layer"  model,  one 
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is  in  fact  dealing  with  infinitely  many  layers  in  the  over- 
burden. Equivalently,  one  may  refer  to  this  as  a randomly 
layered  model. 


Sample  Realizations 

It  is  of  considerable  interest  to  the  interpreter  of 
field  data  to  know  the  extent  to  which  the  variability  in  the 
conductivity  profile  affects  the  apparent  resistivity  and 
kernel  function  curves.  With  this  in  mind,  some  representative 
sample  realizations  of  the  random  conductivity  function  and 
the  corresponding  kernel  and  apparent  resistivity  curves  will 
be  examined. 

Kernel  function:  The  equation  for  the  stochastic  kernel 
in  equation  (2-24)  can  be  rewritten  explicitly  in  terms  of 
z as  follows: 


KjU)  = - 


/ n(z)e 


-2Az 


4k  Ae 


/ n(z)  cosh  2Az  dz 


(4-  1) 


It  is  recalled  from  equation  (2-22)  that  the  deterministic 
two-layer  kernel  is 


Kq(A) 


The  total  kernel  is 


Kq(A)  + eKx(A) 


K(A) 


T-1950 


53 


As  X 


The  asymptotes  of  the  kernel  functions  are  worth  examining. 

* 0 ■ * fct 


and 


as  X ■*  « , K -*•  1 

o 


Obviously,  ■*  0 as  X -*■  0 except  when  the  basement  is 

infinitely  resistive  (k  = 1)  in  which  case  the  expression 
increases  without  limit  and  the  solution  is  not  applicable. 

To  derive  the  other  asymptote  the  equation  for  needs  to  be 
rewritten  as 


KX(X) 


2 (l+ke~2Xh)X 
(l-ke-2Xh) 


h 

/ n(z)e 
o 


-2Xz 


dz 


- 2k  X 

d-ke-Wh,2 


h 

/ 

o 


{e-2X(2h-z)  + 


j~2X(2h+z)}  n(z)  d2 


(4-  2) 


For  large  X (>  5)  the  second  integral  is  negligible  compared 

to  the  first.  (For  X = 4 the  first  integral  is  3 orders  of 

magnitude  larger  than  the  second;  for  X = 5 it  is  4 orders 

of  magnitude  larger.  The  ratio  of  the  first  to  the  second  is 
2X 

approximately  e for  large  X . ) In  addition,  for  large  * , 
the  first  term  in  equation  (4-2)  approaches  -n(0).  Hence, 

K(X)  -*■  as  X -*-0  and  K(\)  -*■  1 -en(0)  as  X -*•  00  . 

The  second  asymptote  is  independent  of  k. 

Some  examples  of  the  sample  realizations  of  the  kernel 
function  K(X)  for  k = ± 0.8  and  c t = 0.05  are  presented  in 
Figures  24  through  27.  Superimposed  on  these  are  the  deter- 
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Figure  25.  Sample  realization  No.  2 of  kernel  function  for 


le  realization  No.  3 of  kernel  function 
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ministic  kernel  functions  Kq(X)  plotted  as  solid  lines  in 
the  figures.  The  random  conductivity  functions  n(z)  are 
shown  in  the  upper  right-hand  corners  of  the  graphs.  The 
parameter  e has  been  set  at  the  relatively  large  value  of 
0.1  to  accentuate  the  deviations  from  the  mean.  The  abscissa 
in  each  graph  is  in  dimensionless  units  of  Xh  . 

These  representative  curves  show  the  range  and  nature  of 
deviations  from  the  mean.  The  most  important  factor  evidently 
is  the  value  n(0).  Depending  on  the  actual  nature  of  n(z) 
the  kernel  functions  may  look  like  those  of  a deterministic 
two-layer  case  (with  appropriate  shifts  in  the  axes)  or,  in 
extreme  cases,  for  example  in  Figure  24,  they  may  resemble  those 
of  three-layer  models. 

The  deviations  in  the  kernel  function  are  more  clearly 
illustrated  if  graphs  of  the  changes  relative  to  the  cor- 
responding deterministic  values  are  plotted.  The  relative 
change  is  defined  as 

eK.(X) 

— = x 100% 

K0(A) 

Curves  showing  these  relative  deviations  comprise  the  upper 
halves  of  Figures  28  to  34,  the  first  four  of  which  correspond 
to  the  same  four  random  functions  used  in  Figures  24  to  27.  The 
asymptotic  values  are  zero  (as  X -*■  0)  and  -tri(O)  x 100%  (as 
X approaches  infinity).  These  deviations  are  either  positive 
or  negative  depending  on  the  actual  form  of  the  random  function 
and  on  the  particular  X values. 
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j Figure  32.  Sample  realization  No.  5 of  percent  deviation  in 

* kernel  and  apparent  resistivity  functions  for 

k - t 0.8. 
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Figure  33.  Sample  realization  No.  6 of  percent  deviation  in 
kernel  and  apparent  resistivity  functions  for 
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Figure  34.  Sample  realization  No.  7 of  percent  deviation  in 
kernel  and  apparent  resistivity  functions  for 
k - ± 0.8. 


T-J.SOU 


66 


It  is  interesting  to  note  that  the  curves  for  k » + 0.8 
and  k = - 0.8  are  identical.  To  see  the  theoretical  basis 
for  this  let  us  write 


I K1 C X ) | h _2  4k2x  h 

illT  " “ n(Z)e  dZ  + 71^*7  o'  nC.).-4“co.h  2Az  dz 

° (4-  3) 

Evidently,  this  ratio  is  independent  of  the  sign  of  k.  The 

magnitude  of  k,  however,  does  affect  the  ratio.  Inasmuch  as 

the  ratio,  for  large  X , depends  only  on  the  first  terra  in 

the  equation  above,  it  is  independent  of  the  magnitude  of  k. 

For  small  X values  however,  the  second  term  becomes  significant; 

its  value  increases  as  |k|  increases.  In  any  case,  this  ratio 

tends  to  zero  as  X approaches  zero.  This  discussion  is 

supported  by  Figures  35  and  36  for  k = - 0.5  and  k = - 0.99 

respectively.  They  are  to  be  compared  to  Figure  28  which  uses 

the  same  realization  of  the  random  function. 

Apparent  resistivity:  As  noted  previously,  the  normalized 

apparent  resistivities  for  the  Wenner  and  Schlumberger  arrays 


are,  respectively, 

00 

Pw*(a)  = 2a  /K( X ) (JQ(Xa)  - JQ(2Xa)}  dX 
0 


(4-  4) 


Ps*(a) 


a2  / K( X ) X Jx(Xa)  dX  . 
0 


(4-  5) 


Only  the  Wenner  array  will  be  used  for  the  following  discussion. 


Since  K( X ) = K (X)  + eK  (X)  equation  (4-4)  can  be  de- 
3 1 

composed  into  deterministic  and  stochastic  parts: 
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w , o 


2a  / KQ(A)  { JQ(Aa)  - Jo(2Aa)}  dA 


(4-  6) 


w . 1 " -4a  / ^z)  / ^0(Xa)  - J0(2Aa)}  dA  dz 


o o 

h 


-8ak  / n(z)  / f2  {J  (Aa)  - J (2Aa)}  dA  dz  (-•- 
o o A ° ° 


where  f1(A,  z) 


A(l-t-ke'2Xh)e~2Az 


(l-ke’2Xh) 


(4-  8) 


f2(X*  2) 


Xe~4Xhcosh  2Az 


(l-ke_2Xh)2 


(4-  9) 


« 

M 

.1 


41 

1 


A rapid  and  reasonably  accurate  procedure  for  evaluating 
Hankel  transform  integrals  is  the  use  of  digital  linear  filters. 
This  method  is  already  well  established  (for  example  Ghosh, 

1971  and  Anderson,  1975).  In  general,  this  approach  is  about 
five  times  faster  than  numerical  integration  by  Gaussian 
quadrature  and  provides  acceptable  accuracy  of  about  four 
significant  figures  (Anderson,  1974). 

In  order  for  this  method  to  be  applicable  however,  the 
functions  to  be  transformed  (fj  and  tz  in  this  case)  must 
be  monotonic  decreasing  as  A * »,  Inspection  of  equation  (4-8) 
reveals  that  f i satisfies  this  requirement  except  when  z * 0. 

In  this  special  case  one  can  expand  fi  as  an  infinite  series, 
find  the  Hankel  transform  analytically  term  by  term  (Erdelyi, 
1954)  and  finally  obtain  a power  series  expression  as  follows: 


i 
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[ fx(X,  0)JQ(Xr)  = 4h  £ 


nk 


q»1  {( 2nh)*  + r2}3/i 


(4-10) 

It  is  interesting  to  note  that  even  when  z • 0 the  digital 
filter  approach  yields  results  of  acceptable  accuracy  when 
compared  with  the  exact  value  as  given  in  equation  (4-10). 

Similarly,  it  can  be  shown  that  ti  approaches  zero  as 
X -*■  ® for  all  values  of  z.  Hence,  we  can  apply  Hankel  trans- 
form filters  for  the  evaluation  of  part  of  equation  (4-7). 

Equation  (4-6)  for  the  deterministic  term  can  also  be 
evaluated  by  the  digital  filter  approach  if  it  is  rewritten 
in  the  form 

-2Xh 


w . o 


1 + 4a  / 


ke 


o 1-ke 


^2Xh 


{ J ( Xa)  - J_(2Xa) } dX  . (4-11) 


Seven  representative  sample  realizations  of  normalized 
Wenner  apparent  resistivity  curves  for  a = 0.05,  e = 0.1  and 
k **  - 0.8  are  shown  in  Figures  37  through  43.  The  abscissae 
are  in  dimensionless  units  of  a/h.  These  curves  are  similar 
to  those  for  the  kernel  functions  (if  the  abscissae  are  re- 
versed) except  for  the  well  known  fact  that  apparent  resistivity 

♦ 

curves  are  not  symmetrical  about  the  line  P_  « 1 for  positive 

W 

and  negative  values  of  k;  the  kernel  function,  on  the  other 
hand,  is  symmetrical  about  the  line  K(X)  * 1 for  positive 
and  negative  values  of  k. 

It  is  intuitively  obvious  that  as  the  electrode  spacing 
a approaches  zero,  the  normalized  apparent  resistivity 
approaches  1 - en(0).  As  the  spacing  a becomes  large  c* 
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Figure  42.  Sample  realization  No.  6 of  apparent  resistivity  function 
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approaches  the  deterministic  value  of  (1  + k)/(l  - k)  as 
the  effect  of  the  random  nature  of  the  overburden  becomes 
negligible.  Most  of  the  information  in  this  part  of  the 
sounding  curve  is  returned  from  the  basement. 

The  relative  deviations  from  the  mean  apparent  resistivity 
curves  are  shown  in  the  lower  halves  of  Figures  28  to  36.  The 
asymmetric  nature  of  the  apparent  resistivity  curves  with 
respect  to  the  sign  of  k is  quite  apparent.  It  is  interesting 
to  note  that  the  curves  for  positive  k are  remarkably  similar 
to  those  for  the  kernel  function. 

As  in  the  case  of  the  kernel  function  an  increase  in 
the  magnitude  of  k increases  the  relative  change  in  the 
apparent  resistivity  curves  for  large  electrode  spacings, 
out  has  no  effect  for  small  spacings.  The  difference  between 
the  curves  for  positive  and  negative  values  of  k is  accentuated 
if  k has  a greater  magnitude. 

Ensemble  Statistics 

It  has  been  well  established  in  the  literature  that  the 
kernel  function  embodies  all  the  information  about  layer  para- 
meters which  are  contained  in  the  apparent  resistivity  curve 
itself.  In  the  interest  of  economy  in  computation  time  the 
study  of  ensemble  statistics  is  confined  to  the  kernel  function 
instead  of  the  apparent  resistivity  function  to  which  it  is 
related  via  a Hankel  transformation. 

Equation  (2-24)  rewritten  in  a more  convenient  form 


becomes 
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KX(X)  = W / (qx  + q2  + q3)  N(u»)  dw 


(4-12) 


where 


W(X) 

_ 4X. 

q1(w. 

X) 

« 2 Xv 

y 

q2  ( k*  > 

X) 

2Xe 

a ^ 

-2Xh 

q3(w, 

X) 

= e 

— 2 X h 1 1 i.2  - 


$ u 


$ y 


The  mean  <Ki(X)>  is  obviously  zero  since  <N(w)>  is  zero.  The 
covariance  is 

2 cd  2 2 

<K1(X1)K1(X2)>  * ^ W(X1)W(X2)  /T(co,X1)T(u),X2)e~a  ^ du) 

/rr  o 

(4-13) 


where  T = qi  + q2  + q3 


Convergence  of  the  integral  is  assured  by  the  presence  of  the 
exponential  term.  The  covariance,  normalized  by  the  standard 
deviations,  is 

<K. (X. )K. (X9)> 

R(X  X2)  - ■ ■■■■-  - J • (4-14) 

{<K1(X1)K1(X1)>  <K1(X2)K1(X2)>} 

The  normalized  covariance,  computed  for  representative 
values  of  X and  a are  shown  in  Figures  44  and  45.  It  is  clear 
that  as  a increases  the  stochastic  kernel  function  becomes 
more  correlated.  In  addition,  the  larger  the  X value  the 
more  correlated  is  the  random  kernel.  Thus,  a more  correlated 
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Figure  44.  Normalized  covariance  of  two-layer  random  kernel 
function  for  a = 0.05. 


LflMB0fl-2 


Figure  45.  Normalized  covariance  of  two-layer  random  kernel 
function  for  a » 0.1. 
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random  function  n(z)  will  be  reflected  in  a more  correlated 
random  kernel  function. 

Sample  Statistics 

Two  hundred  sample  realizations  of  the  random  kernel  as 
given  in  equation  (4-1)  were  computed  and  the  distributions 
tabulated  in  the  form  of  histograms.  Some  of  these  are  shown 
in  Figures  46  and  47.  Using  the  chi-square  goodness  of  fit 
test  all  the  distributions  generated  for  various  a and  X 
values  were  found  to  be  Gaussian  at  the  5%  level  of  significance 

The  computed  means  are  generally  much  less  than  0.05 
while  the  standard  deviations  are  approximately  in  the  range 
0.5  to  1.0.  The  coefficients  of  skewness  are  in  general  less 
than  0.2  in  magnitude,  being  about  evenly  distributed  between 
positive  and  negative  values. 

It  is  known  that  a Gaussian  process  preserves  its  Gaussian 
form  after  a linear  transformation  (Freeman,  1958).  Inasmuch 
as  the  Hankel  transformation  of  the  kernel  function  to  the 
apparent  resistivity  function  is  a linear  process  (Ghosh,  1971) 
one  may  conclude  that  the  random  apparent  resistivities  are, 
at  the  5%  level  of  significance,  normally  distributed  about 
their  mean  values. 

Effect  on  Interpretation  of  Data 

It  is  appropriate  at  this  point  to  pose  the  following 
questions.  What  effect  does  a random  distribution  of  earth 
conductivities  have  on  resistivity  sounding  curves?  To  what 
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extent  is  the  interpretation  of  resistivity  data  affected? 

How  can  these  effects  be  taken  into  consideration  during 
data  interpretation?  These  and  related  questions  will  be 
discussed  in  the  following  paragraphs. 

The  effect  of  the  random  nature  of  the  distribution  of 
conductivities  on  an  apparent  resistivity  sounding  curve  is 
greatesj  for  small  electrode  spacings.  With  small  spacings 
we  are,  effectively,  probing  only  the  shallower  regions  of  the 
earth.  It  is  evident  from  the  results  presented  previously 
that  as  the  electrode  spacing  a decreases  the  apparent 
resistivity  approaches  the  resistivity  of  the  near-surface 
zone,  that  is, 

p(o)  » p°(l  - en(o)  ) 

As  the  spacing  a increases  more  and  more  information  is 
returned  from  the  deeper  regions  of  the  earth  and  the  apparent 
resistivity  approaches  the  resistivity  P2  of  the  lower  layer 
(basement).  In  this  case,  the  random  conductivity  profile 
in  the  overburden  has  a negligible  effect.  Thus,  the  net  effect 
of  the  random  conductivity  distribution  in  the  overburden  is 
to  shift  the  part  of  the  sounding  curve  for  small  a up  or 
down  (in  the  sense  of  the  vertical  apparent  resistivity  axis 
of  the  curve)  relative  to  the  part  of  the  curve  for  large  a 
according  as  p(c)  is  greater  or  less  than  the  mean  overburden 
resistivity  p°.  This  relative  shift  can  be  considerable 
depending  directly  on  how  much  p(c)  deviates  from  p°.  As 
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a practical  example,  the  near-surface  zone  of  a sedimentary 
section  is  probably  drier  than  the  lower  zones.  This  results 
in  an  increase  in  its  resistivity  relative  to  the  rest  of  the 
section. 

The  variability  in  the  resistivity  profile  tends  to  be 
smoothed  out  in  the  sounding  curve.  Consequently,  the  earth 
may  be  viewed  as  a high  frequency  filter  for  the  noise  in 
the  system.  The  shape  of  the  resulting  sounding  curve  depends 
to  a large  extent  on  the  function  n(z).  However,  the  value 
n(o)  has  a disproportionate  effect  on  this.  The  apparent 
resistivity  curve  may  resemble  a two-layer  curve,  as  is  usually 
the  case  (for  example  Figures  39  and  40)  or  there  may  be 
distinct  indications  of  an  "intermediate  layer"  as  in  Figure 
37.  Such  deviations  are  usually  small  and  are  difficult  to 
detect.  In  addition,  measurement  noise  (which  includes 
instrument  noise,  human  error  and  lateral  inhomogeneities) 
tends  to  mask  these  small  deviations.  In  practice  these 
curves  would  tend  to  be  considered  as  two-layer  curves  with 
slightly  different  parameters. 

For  a semi-quantitative  estimate  of  the  magnitude  and 
type  of  errors  introduced  in  the  interpretation  of  such 
curves  we  will  resort  to  standard  curve-matching  techniques. 
Several  two-layer  apparent  resistivity  curves  were  traced 
from  the  standard  curves  of  Orellana  and  Mooney  (1966). 
Deviations  of  -10%  and  -20%  were  then  introduced  into  the 
parts  of  the  curves  which  correspond  to  small  electrode 


